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Introduction 

The purpose of this f port Is to give a brief aunnary of the progress 
■ade In the study of mixed time Integration method, for transient thermal 
analysis of structufs during the period from October 1» 1981 to January IS *1982. 
A simple and Illustrative example problem Is also proposed (as suggested by 
Dr. Olsen) to demonstrate the practicability and usefulness of our proposed 
approach. We believe that «e have now developed an efficient solution 
procedure for predicting the thermal behavior of aerospace vehicle 
structures. Currently, a 2D finite element computer program Incorporating 
these methodologies Is being Implemented. The performance of these mixed time 
finite element algorithms can then be evaluated employing the proposed example 
problem. Related discussion is In later section. 

Sumaary of Developments 

A family of mixed time Integration schemes In which different time 
Integration methods (ln^llclt/expllelt) with different time steps can be used 
In each element group is developed. Its underlying theoretical framework for 
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l«rg« scale aulel-dlBsaslonsl lliwsr cranslsot analysis is studied* The 

effective coeputer iapleaentation aspects of these proposed techniques are 

also investigated* The "active column equation solver" is the key to success 

of this technique* The equations system for each element group is constructed 

in such a way chat they arc uncoupled and hence each group can be integrated 

at Its own group time step with different integration methods* The stability 

characteristic of these mixed time partition algorithms is also studied* An 

energy balance technique is e^iloyed to carry out Che stability analysis* It 

« 

is found chat the critical time step restrictions are governed by each 
individual explicit element group only* Of more significance, these mixed 
time ImpllciC-expllclC methods provide a natural framework for Che further 
development of efficient, clean and modularized computer codes* 

The above developments are summarized in the attached two papers* 
Evaluation of the Mixed Time Implicit^Explicit Algorithms 

A number of simple 10 numerical examples are presented in the two papers 
to demonstrate the feasibilities (i*e* accuracy and stability behavior) of Che 
algorithms* We have formulated herein two finite element thermal models for 
Che Thermal Protection System (TPS) of the space shuttle to evaluate Che 
performance of our mixed time Implicit^explicir. method* The first one is a ID 
model while Che second one is a 20 model* They are depicted in figures 1 and 
2 respectively* The initial phase of the research project is the development 
of the mentioned methodologies in which linear assumption is employed* 
Nonlinear effects such as incemal and external radiation, and natural 
convections are ignored* Further only the assumed mean temperature material 
properties of the varic ts components of the TPS are considered* Its siaterial 
properties are tabulated in Cable 1* The number of elements, Che minimum 
ctiaraccecisclc length I^q, the estimated explicit critical time step ^'mln* 


c f.V, 

^ 1 . ■ ^ 

the pcopoeed laccgraclon oethod end the proposed elenent group clae step for 
eech group ere included in teble 2. Ae cen be seen froa teble 2. due Co the 
various therael ciae scales (e.g. AC . - 0.083 for AL, end 

Utl 

133*6 for RSI), e single incegreclon aechod is definicel 7 not 
effective. For exeaple, if en explicit aethod is eaployed, e tiae step of 
0*083 has to be used; while if en iaplicit aathod is eaployed, there is no 
scebility-iaposed liaitation on At ; however, wide->bandiag and/or non- 
convergence of the resulting aatriz equations aay decrease its advantage. The 
faaily of aixed tiae integration scheaas developed is best suited for this 
type of problea* The attributes of the various tiae integration aethods are 
fully achieved using the proposed approach as can readily be seen froa 
cable 2. It should also be observed that 6C^^ can be sec as high as 33.42 
in this aixed aodel. Subsequently At « 33*42 is eaployed for eleaent groups 
1 and 2 (042 Coating and RSI), At " 66.84 for eleaent group 3 (RSI), 
and At ■ 133.68 for eleaent groups 4 and 5 (RTV, FELT, AL and AIR). However, 
judging froa Che heat loading versus ciae curve (see figure 3) it is advisable 
to pick At . " SI. 4/4 • 12.85 during the heat loading period (i.e. tiae * 0 

lUtI 

to 1200) and At . ■ 33.42 in later tiae. Of course, an auCoaaClc selection 

oiin 

of Che step size based on truncation error and accuracy considerations is 
desirable. It is not being considered here because of the Halted tine frane 
and budget of this investigation. 

The laportance of this proposed nixed tine Inplicic-expllcic finite 
eleaent concept can further be visualized if a nonlinear analysis of the above 
finite eleaent nodels is assuaed. The themal responses of the various 
coaponents of the TPS nay be divided into regions of slowly and rapidly 
varying ceaperatures and they are: 
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1. 042 Coating 

Th« th«nul diffusivlcy Is slaosc ladspsndsat of esapscseurs sod ch« 
ssclascs 4t is so saall (0.25) chat sxpllcic cslculseioa is aoc cose 
sffscClvs sc all. lapllcic calculaCioa Is base suicad sloes only 
llalcad partial caforaaclon and rsfaccorlsaclon of this alsasoC 
group’s affacClva selffoass will ba oaadad. This Is ona of tha aajor 
advaccagas of cha alxad ciaa lapllclt-axpliclc concapc. 

2. RSI 

Tha ehacaal diffusivlcy ehaogas rapidly with caaparaeura» charafora 
raforaaclon and rafaccorlsaclon of this alaoanc group’s affacelva 
sclffnass ara fraquantly raqulrad If lapllcic aachod Is aaployad. 
Hoo^coovarganea and/or wlda-bandlng of Cha rasuleing alaaanc group 
affacc'va sclffnass aay daeraasa cha advancaga of applying cha 
lapllcic MChod hara. Tha aselaacad 4C (basad on aaan Caaparatura 
valua) is 33.42; Charafora axpllclc Is bast suited (no aacrlcas 
oparaclons). 

3. RT7 

Its Charaal diffusivlcy Is Indapandanc of caaparatura and Cha 
asclsMCad &C Is too saall (0.213) chat lapllcic owchod Is 

racoaasndad. 

4. ^ 

Its Charaal diffusivlcy Is fairly Indapandanc of caaparatura (ac 
laasc In cha oparaclonal caaparatura range of Che space shuttle) and 
Its esclaaced 4C Is so stringent (0.083) chat li^llclc aechod Is 
ags'n raconoanded. 

5. FELT 

Its chermal diffusivlcy Is Indapandanc of ceaperacura and Its 
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Mtiwecd &c is Isrgs (51.4), chsrsfors inplicit MChod or sxpllclt 
Mthod csn bs used. For this sxsaplt, we propose Co e^loy iapllclc 
Mthod. However for 3D eslculsCion, in order to reduce the bend** 
width end Che utmecessery aonllaesr cslculstloas (due to Che feet 
Chet the edjecent groups sight be Ispllclt groups too), explicit 
sMChed will be recossended. 

AIR 

The verletioos of the Chersel dlffuslvlty with tesperscure ere saell 
in our rsnge of Interest. Also the estlseted It Is ssell, therefore 
ispllclt sethod is egeln well suited. 

In Che ebove dc celculetlons, I . ■ sin (i } where i end 

sin ^ X y' X y 

ere defined as follows: 

IT 

‘r 



H *,-H| 

2 

The estlseted At Is defined to be A * eonservetlve 

critlcel else step celculeclon. Note Chet s Is the chersel 
dlffuslvlty. 

By virtue of the feet thee each elescne group's effective stiffness 
Is uncoupled to the global essesblled setrlx equetions systes, any 
eleawnt group can be reforsed and refactorized at any instant If 
required without affecting the global equations system. This partial 
factorization procedure can further be enhanced if it is to be 
combined with an iterative update procedure. 
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(4) In th« cn«« of short tla* duration cransiants, oixad tin* explicit- 

explicit is aost effective. It is even more attractive if the nixed 
tine explicit-explicit procedures are unconditionally stable. 

Conclusions end Future Directions 

The work suanarlxed here is not complete. The aost iaportant work in 
progress is the developaent of nixed tlae integration nethods for transient 
thermal analysis of structures suitable for incorporation into aost finite 
eleaent coaputer codes. These alxed tine procedures are currently being 
Integrated into a pilot finite eleaent coaputer program. Detailed studies of 
solutions of the proposed examples are planned. 
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Tibi* !• Neao TMparatur* Material Properties 


MATKKIAL 

042 

COATING 

ESI 

PKOPKRTIES 



X w/e-k 
condiicClvtCy 

1 

I _ 

1.4338 

0.1354 

P kg/«^ 
density 

1665.92 

144.17 

C J/kg-k 
specific heat 

1317.96 

1255.20 

j 

X 1 , 

a- ^ ■ /see 

6.53*10“^ 

7.48x10“^ 

Che real 
dlfCualvlCy 





Ibe nean teaperature aaterlal properties are coeputed from Che average of Chose aC 
12U<) k. 9SU k. 500 k, 477 k, 333 k and 300 k. 
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Iteble 2* Characteriatlc Length and Tine Scales 


MATBKIAL NUMBER OP BSTINATBD INTEGRATION 

TYPB ISLEMENTS L^^„(c«s) At MBTtlOD 


o,^«ENT emit 

TIME STEP 


0A2 

coating 



133.68 


total nunber cf elements •> 240 


1 “ Implicit Integration 
E “ explicit integration 


































